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MASSIVE PLANET MIGRATION: THEORETICAL PREDICTIONS 
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ABSTRACT 

We quantify the utility of large radial velocity surveys for constraining theoretical models of Type II 
migration and protoplanetary disk physics. We describe a theoretical model for the expected radial dis- 
tribution of extrasolar planets that combines an analytic description of migration with an empirically 
calibrated disk model. The disk model includes viscous evolution and mass loss via photoevaporation. 
Comparing the predicted distribution to a uniformly selected subsample of planets from the Lick / 
Keck / AAT planet search programs, we find that a simple model in which planets form in the outer 
disk at a uniform rate, migrate inward according to a standard Type II prescription, and become 
stranded when the gas disk is dispersed, is consistent with the radial distribution of planets for orbital 
radii 0.1 AU < a < 2.5 AU and planet masses M p > 1.65 Mj. Some variant models are disfavored 
by existing data, but the significance is limited (~95%) due to the small sample of planets suitable 
for statistical analysis. We show that the favored model predicts that the planetary mass function 
should be almost independent of orbital radius at distances where migration dominates the massive 
planet population. We also study how the radial distribution of planets depends upon the adopted 
disk model. We find that the distribution can constrain not only changes in the power-law index of 
the disk viscosity, but also sharp jumps in the efficiency of angular momentum transport that might 
occur at small radii. 

Subject headings: solar system: formation — planets and satellites: formation — planetary systems: 
formation — planetary systems: protoplanetary disks — accretion, accretion disks 



1. INTRODUCTION 

Radial velocity and transit searches for extraso- 
lar planets have detected in excess of 170 low- mass 
companions around nearby, mostly Solar-type, stars 
(|Butler et al.l [2006). These detections, which result 
from radial velocity surveys targeting a few xlO 3 
stars, have allowed for an initial determination of 
the distribution of massive extrasolar planets with 
mas s, orbital radius, eccentricity, and stellar metallic 
itv (ISantos. Israelian fc, Mavori [20041: iMarcv et all 1200 
iFischer fc Valentil 120051 : ISantos. Benz fe Mayor! 1200 ' 
The statistical (and, hopefully, systematic) errors on 
these determinations will improve as ongoing surveys 
press to larger orbital radii and smaller planet masses. 
Substantially larger radial velocity surveys of 10 5 to 10 6 
stars, with precision in the ^10 ms _1 range, are tech - 
nically possible over the next decade (|Ge et al.l l2006f ). 
Given this rapid observational progress it is of interest 
to ask how much information - about planet formation, 
planet migration, and the protoplanetary disk - is re- 
tained in the statistical properties of extrasolar planets to 
be potentially tapped via an expansion of existing planet 
samples. Put more bluntly, is it worth obtaining much 
larger samples of planets with properties similar to those 
already known, or does the primary scientific interest for 
future surveys lie in exploring entirely new regimes of 
parameter space? 

At sufficiently small orbital radii, massive extraso- 
lar planets very probably migrated inward from for- 
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mation sites further out rather t han fo rming in situ 
flLin, Bodenheimer fc Richardson] 119961: ITrilling et al.l 
Il998t iBodenheimer. Hubickvi fc Lissauerl [2000D . There 
remains some uncertainty in quantifying 'sufficiently 
small', but it seems likely that massive planet 
formation is most common outside the snow line 
([Hayashi. Nakazawa fc N akagawa 1985h. Protoplanetary 
disk models clearly show that the radius of the snow 
line changes dramati cally with time as the disk evolves 
(|Garaud fc Lin|[2007t ). so to quote a single radius is po- 
tentially misleading. However, for a Solar type star, the 
apparent presence of hydrated minerals in Solar System 
asteroids allows an empirical determination of the lo- 
cation of the snow lin e at a radius of around 2.7 AU 
(jMorbidelli et al.ll2000b . This suggests that most of the 
extrasolar planets currently known, which orbit within 
a few AU of their host stars, derive their properties 
largely via migration. For massive planets, the appro- 
priate regime of migration is thought theoretically to be 
the Type II regime — which inv olves migration within 
a gap in the protoplanetary di sk (jGoldreich fc Tremaind 
1980; Lin & Papal oizoul Il986h — rather than the gap- 
less Type I regime app ropriate to E arth mass planets 
and giant planet cores (|W^ll99l . Since direct ob- 
servational evidence of migration is currently lacking, 
the only tests possible of this theory come from statis- 
tical comparison with the observed properties of extra- 
solar pl anetary systems. Indeed, prior work along these 
lines b y lArmitage et al.l (l2002h , ITrilling. Lunine fc Bend 
(|2002f) and llda fc Lin! (|2004f ) has shown that the distri- 
bution of o bserved planets in orbital radius (and, in the 
case of the llda fc Linl (|2004l ) study, planetary mass) is 
broadly consistent with theoretical expectations based on 
disk migration within an evolving protoplanetary disk. 
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In this paper, we develop more refined predictions for 
the radial distribution of massive planets based on a sim- 
ple analytic model for Type II migration. Our main 
goal is to determine, at least in principle, what might be 
learned from comparisons of large planet samples with 
theoretical models. In §2 we describe the adopted model 
for the protoplanetary disk, and how migration of mas- 
sive planets within the disk is treated. In §3 we compute 
the predicted distribution of planets in orbital radius. 
We compare the pre dictions to the observed distribu- 
tion of planets in the iFischer fc Valentl (|2005f ) sample, 
which has previously been used to study the dependence 
of planet frequency on host metallicity. This sample has 
a clearly specified selection limit, which allows for a reli- 
able statistical comparison between the model and obser- 
vations. In §4 and §5 we investigate the extent to which 
migration leads to a radial variation in the exoplanet 
mass function, and how sensitive it is to structure within 
the protoplanetary disk. These Sections are primarily 
forward-looking, since existing data are too limited to 
support or refute the theoretical model. We conclude in 
§6 with some discussion of the results. 

2. TYPE II MIGRATION WITHIN THE PROTOPLANETARY 

DISK 

We consider a model in which massive planets form 
within an evolving protoplanetary disk at radii be- 
yond the snow line, and migrate inward within a gap 
('Type II' orbital migration). Migration slows and even- 
tually ceases as the gas disk is dispersed. We assume 
that photoevaporation causes disk dispersal. The result- 
ing distribution of massive planets in orbital radius then 
depends upon the disk model (which is reasonably tightly 
constrained by observations); the migration rate (which 
is known reasonably well theoretically); and the rate of 
planet formation in the disk as a function of time. The 
latter can in principle be predicted from a model of mas- 
sive planet formation, but here we treat it as a free func- 
tion. 

The surface density £ of a protoplanetary disk that 
evolves under the combined action of an effective kine- 
matic viscosity v, and mass l oss per unit area E w ind('')) 
is described by (|Pringldll981h . 
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provided that the mass lost in the wind has the same spe- 
cific angular momentum as the disk at the launch point. 
Approximating the angular momentum transport in the 
disk as a time-independent kinematic viscosity that is a 
power-law in radius, 

v = is r*> (2) 

equation ([T]) admits a comp act self-similar solution if the 
mass loss rate is neglig ible (|Lvnden-Bell fc Pringld[l974l : 
lHartmann et al.lll998h . In this solution the surface den- 
sity evolves according to, 

C 



E(fl,T) = 



37W S calci? 7 



T -(5/2-7)/(2-7)exp[-#(2-7) /T ] 



where the scaled variables R and T are defined via 
fiducial radius r sca i e , 
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Fig. 1. — The time dependence of the accretion rate on to the 
central star plotted for different disk models. The three dashed 
curves in the upper panel show the predicted evolution in self- 
similar models with different power-law exponents for the disk vis- 
cosity: 7 = 0.5, 7 = 1.0 and 7 = 1.5. The solid curves show the 
result inclu ding photoeva poration, which is here modeled using the 
solution of Rudcn ( 2004) assuming a disk dispersal time of 6 Myr. 
The analytic cut-off is derived assuming 7 = 1.0, but should be a 
good approximation for general power-law viscosity profiles. The 
lower panel shows the surface density evolution at a disk radius of 
1 AU (the 7 = 1.5 model is the uppermost curve). 
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The constant C is the mass accretion rate at t = as 
r — > 0. Other quantities, such as the mass accretion 
rate M and radial velocity v r in the disk, can be derived 
straightforwardly using these expressions. 

We constrain the parameters of the similarity solution 
(equation[3]) to be consistent with the observational anal- 
ysis of the ti me dependence of T T auri accretion rates 
presented by iMuzerolle et alj |2000). We fix all mod- 
els to have an initial accretion rate M(t — 0) = 3 x 
1CT 7 Afoyr" 1 , which decays to M = 3 x 1CT 10 Moyr" 1 
in 10 7 yr if the disk lives that long before being dis- 
persed. We also fix r sca i c = 10 AU. Table 1 shows the 
values of viscosity normalization constant uq, and initial 
disk mass Mdi s k, that meet these constraints for differ- 
ent values of power-law index 7. Not surprisingly, since 
the basic observational constraints have been known for 
some time, the disk models we favor have parameters 
similar to those considered in many previous studies 
(lHartmann et all 119981: I Armitage. Clarke fc Pallal 120031 : 
lAlexander. Clarke fc Pringld I2006D . Figure [1] shows the 
resulting evolution of the accretion rate for these models. 

The power-law decline in the late-time accretion rate 
implied by the similarity solution does not yield the sharp 
transition between accreting Classical T Tauri stars and 
non-accreting Weak-L i ned T Tauri stars tha t is observed 
(|Simon fc Pratol[T995t IWolk fc Waited [l99l . It is plau- 
sible that this transitio n is driven by photoevaporation 
(|Ballv fc Scovilldfl982l) fr om the outer disk, which acts 
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Power law index of viscosity 7 Viscosity normalization vq Initial disk mass 

05 1.52 x 10 s cm 3 / 2 s- 1 0.057 M Q 

1.0 15.62 cms" 1 0.051 M Q 

1.5 1.58 x 10" 6 cm 1 / 2 s- 1 0.058 Mm 



TABLE 1 

Properties of self-similar disk models that meet the observational constraints 
on the evolution of the stellar accretion rate discussed in the text 



to starve the inner disk (thereby allowing it to drain vis- 
cously onto the star on a short time scale) once the accre- 
tion rate becomes comparable to the wind mass loss rate 
(|Clarke. Gendrin fe Sotomavorl l2001). If photoevapora- 
tion is driven by irradiation from the central star, then 
the simplest analytic model s predict that the mass loss 
rate per unit area scales as ()Ho llcnbach et al1ll994f ). 

Swind oc r- 5 ' 2 (8) 

exterior to some critical radius r; n , with zero mass loss 
from smaller radii. The critical radius is given to order 
of magnitude by the radius where the sound speed in 
photoionized gas (T ~ 10 4 ) first exceeds the local escape 
speed. This is a few AU for a Solar mass star. The nor- 
malization of the mass loss rate depends upon the square 
root of the ionizing flux $, which is hard to measure ac- 
curat ely but whose value can be constrained observation- 
ally ([Alexander. Clarke fc Pringld l2005t iPascucci et al.l 
12001 . Much more detailed numerical models of 
photoevaporation are now availab le (|Font et alj l2004t 
lAlexander. Clarke fc Pringld l2006f) , but the additional 
complexity they involve is not warranted for this ap- 
plication. Accordingly, we assume a time-independent 
mass loss rate of the form given by equation ([8]), with 
ri n = 5 AU, and zero mass loss at smaller radii. 

Once photoevaporation is included (in this simplified 
form), it is still possib le to derive a Green's function so- 
lution to equation ([1]) (|Rudenll2004D . but there is no com- 
pact form for the evolution of £(r, t) analogous to equa- 
tion ©. At late times (i.e. after the disk has evolved for 
several viscous times) and small radii (r < r; n ), however, 
it is possible to derive an expression for the reduction in 
the inner accretion rate and surface density caused by the 
wind. Defining to a s the time at which the inner accre- 
tion rate falls to zero as a consequ ence of the mass loss 
from the outer disk, iRudenl (|2004l ) finds that the time 
dependence of the accretion rate at small radii follows, 



M = 



3/2 



self — similar? 



(9) 



where Mseif— similar is the accretion rate evolution pre- 
dicted by the self-similar model in the absence of any 
mass loss. The suppression of accretion implied by 
the term in parenthesis is derived for 7 = 1, but 
should also be approximately valid for other values of 
the power law viscosity index. Using numerica l solu- 
tions to equation , we have verified that the IRudenl 
(|2004j ) formula provides a good description of the evo- 
lution of the accretion rate and inner surface density 
during the transition, and hence we apply the same cut- 
off to generate the surface density evolution including 
wind loss E(r, f) from the self-similar prediction (equa- 
tion [5]). Base d on observations of disk frequency in 
young clusters (|Haisch. Lada fc L ada 2001), we take to = 



6 Myr, which yields an accretion rate evolution shown 
as the solid curve s in Figure [J As has been empha- 
sized previously f Clarkc, Gcnd rlnfc Sotomavorl 120011 : 
lArmitage. Clarke fc Pallal 12003) . to a good approxima- 
tion the disk evolution proceeds as if there is no mass 
loss while M 3> Mwindj and is then dispersed rapidly 
once mass loss becomes significant. Figure [T] also shows 
the evolution of the inner (1 AU) disk surface density. 
Irrespective of the value of 7, all of the disk models con- 
sidered have an initial surface density that is substan- 
tially in ex cess of the minimum m ass Solar Nebula refer- 
ence value (|Weidenschill iug 1977), but this drops rapidly 
as the disk evolves. By the time that massive planets 
migrate through this region (after several Myr of evolu- 
tion) the surface density has fallen by at least an order of 
magnitude. The import of this is that the local disk mass 
during migration is typically smaller than the mass of gi- 
ant planets, a regime which leads to significantly slower 
inspiral. 

Sufficiently massive planets are able to open a gap 
within the gaseous protoplanetary disk and, having 
opened a gap, migrate i nward 3 in lockstep with the gas 
via Type II migration (|Lin fe Papaloizoul 119861 ). For a 
planet at orbital radius r p , the gravitational torques ex- 
erted by the planet on the surrounding disk are strong 
enough to maintain a gap provided that the mass r atio 
q = Mp/M* satisfies (Tak euchi. Mivama fc Linl fT996). 



> 



1/2 

V f 
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where c s and Q p are the sound speed and angular ve- 
locity in the disk at the radius of the planet, and a is 
the dimensionless Shakura -Sunyaev viscosity parameter 
(Shak ura, fc Sunvaevlll973T ). Noting that the thickness of 
the disk h ~ c s /f2 p , we find that the 7 = 1 disk model 
specified in Table 1 would correspond to an equivalent 
a ~ 5 x 10~ 3 (matching at 5 AU around a Solar mass 
star assuming that the disk has h/r = 0.05). The viscous 
gap opening criteria is then satisfied for planets of mass 
exceeding about 0.2 Mj. The ther mal gap opening con- 
dition is satisfied at a similar mass (|Bate et al.ll2"00l . In 
this paper, we exclusively consider planets with masses 
of Mj and above, which should accordingly be safely in 
the gap-opening regime of parameter space. 

To compute the Type II migration rate within the 
disk model specified above, we draw on the results of 
ISver fc Clarkd (|1995l ). These authors noted that the rate 

3 Here, we consider exclusively planets that form close 
enough in that the sense of gas flow and migration is in- 
ward — planets that form further out may instead absorb 
the a ngular momentum of t he inner dis k and move out- 
ward l|Veras fc Armitagd 12004 : IMartin et~aTI 120071 ). or become 
stranded outside the annula r hole that is predicted to form as 
photoevaporation proceeds ( Clarke. Gendrin & Sotomayoj] 120011 : 
Matsuyama, Johnstone & Murray 2003). 
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Fig. 2. — The critical accretion rate above which a 1 Mj planet 
migrates as a test particle is plotted as a function of radius for the 
different disk models. 

of migration, which is sometimes assumed to equal the 
radial velocity that gas in the disk would have in the 
absence of a planet, is suppressed once the planet mass 
exceeds a local estimate of the disk mass. Specifically, 
they defined a measure of the importance of the disk 
relative to the planet, 



47r£r2 
B = ^ 



(11) 



which is small if the planet dominates the angular mo- 
mentum budget of the planet + disk system, and large 
otherwise. The Type II migration rate is then given for 
B > 1 (the disk dominated case) by, 

f p = v r (12) 

while for B < 1 (the planet dominated limit), 

r p = B 1/2 v r (13) 

where we have assumed (consistent with equation [2]) that 
the efficiency of angular momentum transport is indepen- 
dent of the surface density. For our disk models, Figured] 
shows the critical accretion rate below which B < 1 for a 
Jupiter mass planet. In all these models, slowdown of the 
migration rate occurs due to the inertia of the planet for 
accretion rates substantially greater than those at which 
photocvaporation starts to influence the disk evolution. 
The effect on the nominal migration time scale, 



tn 



(14) 



is shown in Figure [31 The B < 1 limit is applicable at 
all radii of interest in a disk with M = 10~ 10 M©yr _1 , 
and within radii less than a few AU (depending on the 
planet mass) in a disk with M = 10 -8 M Q yr _1 . 

Once planets exceed the critical gap opening 
mass by a significant factor, the rate of accre- 
tion onto the planet across the g ap drops rapidly 
(|Lubow. Seibert &: Artvmowicd Il999f) . Our goal is to 
compare theoretical predictions for the planet distribu- 
tion with a complete subsample of extrasolar planets, 



Fig. 3. — The predicted migration time scale as a function of 
radius in the 7 = 1 disk model. The solid lines show the migration 
time scale for planets of masses 0.5 Mj, 1 Mj, 2 Mj and 4 Mj 
when the disk accretion rate is 10 — 8 MQyr - 1 (more massive plan- 
ets have longer migration time scales at small orbital radii). The 
dashed lines show results for the same mass planets in a disk with 
M = IO" 10 Moyr" 1 . 

which requires a cut at approximately 1.5 Mj. Most 
of the planets that survive this cut have masses substan- 
tially greater than the gap-opening mass, so as a first ap- 
proximation it seems reasonable to ignore the possibility 
of mass accretion during Type II migration, and assume 
that negligible mass is accreted across the range of radii 
(interior to 2.5 AU) over which w e make the comparison. 
The results of iBate et al.l ((2003), which show that plan- 
ets can grow rapidly to masses beyond that of Jupiter 
while suffering little migration, support this approxima- 
tion. It is then easy, using equations (fT2")) and (|13[) . to 
calculate the final (following disk dispersal) orbital radius 
of a planet that forms with mass M p at radius Of orm and 
time £form in one of the disk models specified in Table 1 . 
The resulting mapping between the formation conditions 
and the final state is the basic input needed to make a 
prediction of the resulting planet distribution. 

The above model represents our attempt to define a 
'best-guess' description of Type II migration within the 
protoplanetary disk. One may note that in the planet- 
dominated B < 1 limit, the predicted rate of Type II 
migration is reduced by the square-root of B rather than 
the full ratio of the local disk mass to the planet mass. 
We will refer to this as partial suppression of the migra- 
tion rate. The rate is not fully suppressed because gas 
accumulates close to the tidal barrier as the disk / planet 
system evolves, increasing the torque be yond the valu e 
that would occur in an unperturbed disk (Pringlc 1991). 
In practice, however, some gas may overflow the gap to 
be either accreted by the pl anet or to flow into an in - 
ner disk interior to the gap (jLubow fc D'Angelol 1200(f ). 
Mindful of this, we have also computed models in which 
the migration rate in the planet dominated regime is fully 
suppressed, so that, 

r p = Bv r . (15) 

The difference between the partially and fully suppressed 
calculations provides an indication as to how sensitive the 
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Fig. 4. — Predicted final radii for 1 Mj planets as a function 
of the time at which they formed at 5 AU in the protoplanetary 
disk. The three closely spaced curves on the right-hand-side of 
the Figure show results for disks with 7 = 0.5 (long-dashed lines), 
7 = 1 (solid lines) and 7 = 1.5 (short-dashed lines) respectively, 
ass uming that migration is suppressed at small radii according to 
the lSverfc "Clar ke (1995) model. The left-hand curves show results 
if, instead, migration is fully suppressed by the ratio of the local 
disk mass to the planet mass. Differences between migration treat- 
ments are evidently more significant that differences between the 
disk models. For the favored models in which migration is partially 
suppressed, the survival time of planets in the disk is in the range 
of 10-20% of the disk lifetime. 

results are to uncertainties in the treatment of massive 
planet migration. 

3. PREDICTED DISTRIBUTION OF PLANETS IN ORBITAL 

RADIUS 

Using the analytic disk evolution and migration mod- 
els described in §2, we compute the radius afi na i at which 
massive planets become stranded as a function of the 
time tform at which they form within the protoplanetary 
disk. The parameters of the model are the planet mass, 
the power-law index of the disk viscosity, the formation 
radius, and whether the migration rate is partially or 
fully suppressed. Illustrative results for a formation ra- 
dius of 5 AU and a planet mass of 1 Mj are shown in 
Figured! As is well known, planets that form at earlier 
epochs migrate to smaller final radii, and as the orbital 
radii decr ease the window of allowed formation times also 
narrows (jTrilling. Lunine fc Ben3 l2002t lArmitage et~all 
2002). For the favored model, in which migration is only 
partially suppressed, all surviving planets must form dur- 
ing the last 1-1.5 Myr of the disk lifetime. This is not so 
short a window as to imply that the planet formation pro- 
cess must be worryingly lossy, but it does imply that the 
final masses of giant plan ets might well reflect details of 
the d isk dispersal process (|Shu. Johnstone fc~ Hollcnbach 
1993). Planets can form and survive across a larger frac- 
tion of the disk lifetime if instead migration is fully sup- 
pressed. There are also significant differences in the out- 
come that depend upon the adopted disk model. These 
arise because of the differing profiles of the radial veloc- 
ity as a function of radius, but they are less significant 
than differences in the treatment of migration. 

To translate these results into a prediction for the or- 



FlG. 5. — The predicted distribution of massive planets as a func- 
tion of radius. All the curves show results for Jupiter mass planets 
(note that the mass dependence is sufficiently weak that essentially 
identical results apply also for higher masses). The solid curve 
shows the results for the favored model, in which planets form at 
5 AU at a constant rate in a 7 = 1 disk, and migration is partially 
suppressed at small radii. The two dotted curves show the effect of 
varying the disk viscosity to 7 = 0.5 or 7 = 1.5. The long-dashed 
curves show the effect of assuming that planets form at 3 AU (the 
curve furthest to the left at 2.5 AU) or 10 AU rather than 5 AU. 
The short-dashed curve is the predicted distribution in the case in 
which migration is fully suppressed. The data points (shown with 
approximate errors) show the number of observed extrasolar plan- 
ets with 1.65Mj < M p sin(i) < lO Mj in a uniform and com plete 
sample of planets constructed from I Fischer fc Valcnti (200jj|). 

bital distribution of massive planets, we note that, 

dN p dt folm dN p 

~a = 1 x j — (16) 

da dafinai dif orm 

where dJVp/dif orm is the rate at which planets of a given 
mass form in the outer disk. The simplest assumption 
is that this rate is a constant, or at least can be ap- 
proximated as such, over the interval of time near the 
end of the disk lifetime during which planets can form 
and survive migration without being swept into the star. 
Making this assumption, we compute and plot in Fig- 
ure [5] diVp/d log a for a variety of models, which differ in 
the assumed formation radius, disk model, and migration 
treatment. Numerical values for three of the these mod- 
els, which assume af or m = 5AU, partial suppression of 
the migration velocity, and 7 = 0.5, 7 = 1 or 7 = 1.5, are 
tabulated in Table 2. These three cases roughly bracket 
the range of predicted outcomes for all of the models that 
we have considered. 
Although all of the curves have the same general form 
- relatively few planets are predicted to be stranded 
at very small orbital radii, with the number increasingly 
rapidly as a increases — there are significant differences 
which may in principle leave signatures in the observed 
distribution of extrasolar planets. Most importantly, the 
predicted number of planets per logarithmic interval in 
semi-major axis rises sharply as the formation radius is 
approached. A model in which planets are assumed to 
form typically at 10 AU rather than 5 AU results in 
a much smaller predicted number of planets at 4 AU, 
whereas assuming a typical formation radius of 3 AU (i.e. 
immediately outside the snow line) leads to a predicted 
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Iog(a/AU) diV(7 


= 0.5)/dlog(a) d7V(7 


= 1 0) /d \oe(a) dN(~< 


= 1 Sydloela) 


-1.0 


0.152 


0.263 


0.420 


-0.9 


0.184 


0.298 


0.423 


-0.8 


0.221 


0.338 


0.489 


-0.7 


0.266 


0.384 


0.529 


-0.6 


0.320 


0.436 


0.573 


-0.5 


0.385 


0.497 


0.623 


-0.4 


0.463 


0.567 


0.679 


-0.3 


0.559 


0.650 


0.743 


-0.2 


0.676 


0.747 


0.816 


-0.1 


0.820 


0.863 


0.899 


0.0 


1.000 


1.000 


1.000 


0.1 


1.226 


1.170 


1.117 


0.2 


1.518 


1.380 


1.261 


0.3 


1.912 


1.656 


1.448 


0.4 


2.458 


2.035 


1.705 


0.5 


3.283 


2.592 


2.080 



TABLE 2 

The predicted number of planets per logarithmic interval in semi-major axis for models in which planet formation occurs 
at 5 AU and the migration rate is partially suppressed. Results for three different values of 7 are quoted. These results 
approximately bracket the shallowest to steepest curves obtained for the models shown in figure [5] results are quoted 
for 1 Jupiter mass planets, but the mass dependence is weak - much weaker than the differences between the different 7 

MODELS. The numbers have been normalized to unity at 1 AU. 



pile-up of planets as that radius is approached. Observa- 
tional detection of a rapid rise in the number of planets at 
a particular orbital radius would then be a signature of a 
preferred orbital radius for giant planet formation. Also 
easily detectable are the differences between the partially 
and fully suppressed migration models. If migration is 
fully suppressed, then the number of planets (expressed 
as diVp/dloga) at 0.1 AU is predicted to be around half 
the number at 1 AU, whereas for partial suppression this 
ratio is around 0.25. This result is easily understood - 
differences in the treatment of migration are most im- 
portant at those small radii, interior to 1 AU, where the 
disk mass is lowest. 

3.1. Sensitivity to the disk dispersal mechanism 

Our parameterization of disk evolution sweeps much 
that is poorly known about photoevaporation into a sin- 
gle parameter - the disk dispersal time to- Physically, 
rapid dispersal (small values of to for a given set of disk 
initial conditions) will occur for large values of the ioniz- 
ing flux 4>, which will drive a stronger photoevaporative 
flow. Quantit atively, the total mass loss rate from the 
disk scales as (Hollcnb ach et al"1ll994f) . 

/ $ \ 1/2 / M \ 1/2 
M wmd . 1.5X10-^) M yr- 

(17) 

where the prefactor has been adjusted from the 
analytic value to better match the numerical re- 
sults of iFont et alj (|2004f) . and the fiducial ionizing 
flux has been cho sen to be consistent with obser- 
vational estimates ((Alexander. Clarke &: Prin gle 2005; 
iPascucci et al.|[2007h . which are, however, subject to sub- 
stantial uncertainties. 

With our disk models, a wind mass loss rate in the 
range between 10 -10 M0yr _1 and 10 -9 M Q yr -1 yields 
the dispersal time of 6 Myr that is our baseline assump- 
tion. Our fiducial parameters are therefore consistent 
with standard photoevaporation models. A consequence 
of this, however, is that the disk mass at the epoch when 
surviving planets form is quite small. Figure [6] shows the 



gas disk mass at the moment planets form as a function 
of the final planet semi-major axis. For the models we 
have considered, the typical disk masses are only a few 
Jupiter masses if migration is partially suppressed (in 
the fully suppressed case, on the other hand, surviving 
planets form substantially earlier when there is plenty 
of disk gas remaining). These small masses directly re- 
flect the efficiency of the migration process - relatively 
modest amounts of gas are able to drive substantial mi- 
gration. However, they do pose a possible consistency 
problem: planets that form in very low mass disks evi- 
dently perturb the disk structure substantially, and this 
might affect the final planetary distribution. 

To explore this possibility, we considered two options. 
First, within the formalism developed here, we have in- 
vestigated the effect of reducing t so that surviving plan- 
ets form earlier when the disk mass is higher. As shown 
in Figure [6j a model computed with to = 3 Myr yields 
an almost identical planet distribution. Second, we have 
compared the r esults for 7 = 1.5 with t he numerical mod- 
els presented in I Armitage et all ((2002) . In the numerical 
models, planet masses are restricted to be no larger than 
a local estimate of the disk mass, and the feedback of a 
massive planet on the structure of a very low mass disk 
is explicitly followed. Reasonably good agreement is ob- 
tained between the numerical and analytic schemes, with 
the magnitude of the differences being comparable to the 
differences between the curves plotted in Figure [5] 

A related question is how important the nature of the 
disk dispersal mechanism is for the resulting planet dis- 
tribution. Physically, a model in which there is no mass 
loss from the disk exterior to the planet (either via a 
wind, or via accretion on to the planet) fails to yield 
any surviving planets - ultimately the disk, no matter 
how small, absorbs the orbital angular momentum of the 
planet and drives it to small radius 4 . However, many 

4 Strictly, this is only true if the inner disk is subject to a zero- 
torque boundary condition. A steady-state solution is possible if 
there is instead angular momentum injection at small radii, for 
example due to the interaction between the disk and the stellar 
magnetosphcrc (Armitage & Clarke 1996). 
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Fig. 6. — The predicted distribution of massive planets as a 
function of radius (upper panel) is plotted together with the disk 
mass evaluated at the formation epoch (lower panel). The solid 
curve shows the results for the favored model, in which planets 
form at 5 AU at a constant rate in a 7 = 1 disk, and migration 
is partially suppressed at small radii. The short-dashed curve - 
almost indistinguishable in the upper panel - shows the effect of 
reducing the disk dispersal time tp from 6 Myr to 3 Myr. The 
long-dashed curve shows the results of assuming an instantaneous 
dispersal of the disk at 6 Myr. 

other models in which the disk is dispersed rapidly do 
yield sensible distributions, so it is of interest to assess 
how sensitive the planet distribution is to the specifics 
of disk dispersal. To gauge this, we have considered 
an extreme model in which the disk evolves viscously 
for 6 Myr without any mass loss, and is then instanta- 
neously dispersed. Figure [5] depicts the resulting planet 
distribution. At small radii (within about 1.5 AU) this 
model tracks the fiducial photoevaporative case closely, 
but further out the instantaneous dispersal model yields 
a much flatter distribution. This reflects the fact that 
in the case of photoevaporation, migration of the last 
planets to form is limited by the rapid drop in the disk 
surface density, and these planets pile up at radii rela- 
tively close to their formation sites. We conclude that 
the innermost part of the extrasolar planet distribution 
(within roughly an AU) ought to be largely independent 
of how the disk is dispersed, but that the distribution 
further out does depend on whether photoevaporation 
or some other mechanism is at work. 

3.2. Effect of the dispersion in disk properties 

The disk models used in this paper have been adjusted 
to approximately reproduce the mean lifetimes and ac- 
cretion rates of observed disks (|Haisch. Lada k Ladal 
l200HlHartmann et al.|[l998h . This procedure ignores the 
fact that, observationally, there is a large dispersion in 
the accretion rate at a given age. This dispersion ex- 
ceeds that expected from measurement uncertainties in 
accretion rates and stellar ages. 

The origin of the intrinsic dispersion in disk properties 
as a function of system age is not known, and hence it 
is impossible to make a blanket statement as to whether 
consideration of a mean model is reasonable or not. Wc 
can, ho wever, distinguish some possibi lities. One pos- 
sibility (|Armitage, Clarke k Pallall2003D is that the dis- 



persion in disk properties arises from a dispersion in disk 
initial conditions (disk mass, disk angular momentum). 
Since the planet distribution depends only upon the disk 
evolution close to the transition epoch (independent of 
the absolute timing), in this scenario the ensemble distri- 
bution averaged over the population would be expected 
to be the same as the mean model we have computed. 
Of course, giant planet formation would be less proba- 
ble in short-lived di sks with low initial surface densities 
(|Pollack et all [19961. so in practice most planets would 
form in those disks with higher initial surface densities 
at radii of a few AU. 

A second possibility is that the observed dispersion in 
accretion rates arises from changes in individual accre- 
tion rates on timescales intermediate between the disk 
lifetime and observable timesc ales. Apart from thermal 
instabilities (|Bell k Linl 11994 ). which are only likely to 
be relevant for initial accretion rates higher than those 
we have considered here, no obvious candidate instabili- 
ties that might yield large fluctuations in M are known, 
but it remains possible that M(r) might vary rapidly. If 
such variability extended across the radial range consid- 
ered here, it is likely that the migration history of planets 
would be altered. 

3.3. Comparison with a uniformly selected dataset 

We compare the models to a subsample of extrasolar 
planets detected by the Keck / Lick / AAT planet search 
program. Initially, we proceed rather conservatively, and 
define a subsample that is as complete and unbiased as 
possible over specified intervals in planet mass and or- 
bital radii. Our procedure is as follows: 

1. We start with a sample of 850 stars of FGK spectral 
types, targeted by the Lick / Keck / AAT planet 
search program, for which 10 or more radial veloc- 
ity measurements over a period of at least 4 years 
are available. Th i s sam ple is listed in Table 3 of 
iFischer k Valentil (|2005l ). For these stars, hypo- 
thetical planets that yield a radial velocity semi- 
amplitude K > 30 ms _1 have nearly uniform de- 
tectability provided that the orbital period is less 
than 4 years. The orbital period restriction corre- 
sponds to a maximum semi- major axis of 2.5 AU, 
assuming Solar mass h osts. From this sample, 
IFischer k Valentil (|2005f) list 46 stars which host 
known extrasolar planets. In some cases, these are 
multiple systems. 

2. In many cases, additional radi al velocity data is 
availa ble beyond that used by IFischer k Valentil 
(2005). We therefore updat e the orbital elements 
quoted bv lFische~ k Valenti (2005) to match those 



reported by iButler et all (|2006l ). Typically the 



changes to the derived masses and semi-major axes 
are rather modest (here we ignore eccentricity, 
which in some cases is more substantially altered). 
Including the multiple planets in some systems, 
this yields a sample of 59 planets. 

3. Finally, we define a subsample that includes only 
planets that are massive enough to be detectable 
across the entire range of orbital radii for which 
the survey is complete. At 2.5 AU, K = 30 ms" 1 
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Fig. 7. — Comparison of the predicted planet distribution with 
the observed one. The solid histogram shows the cumulative radial 
distribution of massive planets with M v sin ( i) > 1 .65Mj. The sam- 
ple is based on that of Fischer & Valcnti (2005), with the planet 
properties updated using the compilation of IButler et al.l (2006). 
The solid curve is the prediction of the fiducial model for planets 
of mass 3Mj formed at 5 AU in a disk with 7 = 1 and partial 
suppression of the migration velocity. The dashed curve shows the 
predictions of a model in which 7 = 1.5 but all other parame- 
ters remain fixed. This, and other alternate models which predict 
flatter distributions, are not supported by the data, though the 
distributions cannot be distinguished at much better than the 95% 
confidence level. 

corresponds to a planet mass of 1.65 Mj, so we dis- 
card planets with smaller masses. We also cut the 
sample at the high mass end, somewhat arbitrar- 
ily, at 10Mj. In terms of radial extent, we include 
those planets with 0.1 AU < a < 2.5 AU. This 
excludes planets with very tight orbits whose dy- 
namics may have been influenced by tidal effects 
and / or penetration into the protostellar magne- 
tosphere. 

It is somewhat striking — given the large number of ex- 
trasolar planets now known — how little of the publicly 
available data can be used for the sort of statistical com- 
parison we are attempting here. Demanding both com- 
pleteness (which restricts the stellar sample used), and 
lack of mass bias (which restricts the minimum mass), 
we are left with only 22 massive planets. Obviously this 
small sample size limits the statistical power to constrain 
theoretical models. However, the sample should be free 
of any systematic biases, which may not have been true 
for earlier analyses that typically used larger compila- 
tions of detected planets. For illustrative purposes, we 
divide the planets in our subsample into four logarithmic 
bins in semi-major axis (between 0.1 AU and 2.5 AU) 
and plot the number of planets in each bin over the the- 
oretical curves in Figure The observed distribution, 
in agreement with the theoretical models, rises rapidly 
with increasing semi-major axis. 

Figure [7] shows the cumulative distribution of orbital 
radii for massive planets in the sample, together with 
selected theoretical curves derived from the differential 
distributions plotted in Figure [5] The baseline model, 
in which planets form at a constant rate at relatively 
small orbital radii (5 AU) in a protoplanetary disk with 



7 = 1, is shown as the solid line. This model is consis- 
tent with the available data (Kolmogorov-Smirnov test 
P = 0.3). Of the variant models considered, those in 
which migration is fully suppressed, 7 > 1, and / or 
planets are considered to typically form at larger radii 
are disfavored, as these models predict relatively more 
planets at sub-AU orbital radii that are not observed. 
A model in which 7 = 1.5, with the other parame- 
ters remaining unchanged, is shown as the upper dashed 
line. This model is disfavored at modest statistical sig- 
nificance (Kolmogorov-Smirnov test P — 0.06). Several 
variant models that we have considered are disfavored 
at roughly comparable confidence levels (95%), but with 
such a small observed sample few definitive conclusions 
can be drawn. 

As noted above, our sample selection procedure is de- 
liberately conservative. We have also compared the the- 
oretical model to a larger sample, obtain ed by taking 
the en tire catalog of extrasolar planets from IButler et all 
(2006) and retaining those planets with orbital radii 
0.1 AU < a < 2.5 AU and masses IMMj < M p sm(i) < 
10Mj (i.e. the same mass and radius cuts as before, 
but relaxing the condition that all planets were detected 
from one survey with well-known selection properties). 
The resulting sample of 53 planets has a radial distri- 
bution that, by eye, looks very similar to the 22 planet 
sample, and the greater numbers give more power to dis- 
criminate between models. The baseline model (7 = 1, 
formation at 5 AU, partially suppressed migration) re- 
mains consistent with the data, while the second model 
discussed above in which 7 = 1.5 is more clearly incon- 
sistent (P = 3 x 10~ 3 ). We do not place much weight on 
th ese results, since th e compilation of planets reported 
in IButler et al.l (|2006h is not uniformly selected. How- 
ever, they illustrate clearly that only modest expansion 
of the sample size — by around a factor of two — would 
allow us (within the formal confines of the model) to test 
whether the observed distribution of massive extrasolar 
planets is or is not consistent with theoretical models 
that differ in their assumptions as to planet formation 
radius, migration rate, and disk properties. Since some 
of these parameters are partially degenerate, somewhat 
larger unbiased samples would be needed to try and pin 
down empirically a single favored model. 

4. RADIAL VARIATION OF THE MASS FUNCTION 

At small orbital radii where B < 1, the migration ve- 
locity is dependent on the planet mass. In the case in 
which migration is partially suppressed, f p oc M p 1 , 
and hence the fraction of time over which the most mas- 
sive planets (10M p ) can form and survive without mi- 
grating in to the star is a factor of a few greater than 
for Jupiter mass planets. In a steady-state disk, how- 
ever, the fractional suppression of the migration velocity 
as a function of radius is fixed for each planet mass (at 
radii where B < 1). Hence, it is approximately true that 
the relative number of planets left stranded at different 
radii is independent of the planet mass, and the predicted 
mass function of gian t planets is independent of radius 
(| Armitage et al.ll2002fl . 

In more detail, however, migration does lead to some 
fractionation of the planet mass function. More massive 
planets must form earlier than their less massive coun- 
terparts if they are to end up at the same orbital radii 
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Fig. 8. — The magnitude of the predicted variation in the ex- 
trasolar planet mass function with radius in different models. The 
solid curve shows results from the fiducial model in which planets 
form at 5 AU, 7 = 1, and the migration velocity is partially sup- 
pressed. The short-dashed curve shows results for 7 = 1.5, with the 
remaining model parameters unchanged from their fiducial values. 
The long dashed curve shows the effect of assuming that a change 
of migration regime occurs between high and low mass planets. In 
all cases the slope of mass function at small radii has been fixed at 
a = 1. 

after disk dispersal, and the evolution of the disk is itself 
time-dependent (being dominated by viscous evolution 
early on, with a change to a steeper decline in the mass 
accretion rate later on due to the influence of photoevap- 
oration). To quantify the extent to which this results in 
a radial variation of the mass function we assume that at 
small orbital radii the mass function can be written as a 
power-law, 



and compute the predicted radial variation in a. Fig- 
ure [8] shows the results for two of the models that we 
have discussed earlier, including the baseline model. For 
dcfiniteness, we set a = 1 at small radii, though this 
choice is arbitrary. The predicted variation in the mass 
function with orbital radius in these models is non-zero, 
but in these and other similar models we have computed 
in which the migration properties of low and high mass 
planets are the same the magnitude of the change is small 
- of the order of Aa = 0.1 or less. It is straight- 
forward to show that variations of this magnitude are 
undetectable with feasible surveys. The conclusion is 
that any detectable change in the mass function with ra- 
dius (at small orbital radii where migration rather than 
mass growth is the dominant effect) would have to be at- 
tributed to effects other than mass-dependent migration. 

There is one exception to this rule. If the migration 
regime for high mass planets is different from that for low 
mass planets, then large changes in the mass function 
with radius result. Such a change in migration regime 
is conceivable — for example it is possible that high 
mass planets completely prevent mass flow across the 
;ap, while low mass planets allow sig nificant mass flow 
Lubow. Seibert &: Artvmowicz| [i999). In this case, one 
would expect less perturbation to the disk structure im- 



FlG. 9. — The cumulative mass function of known extrasolar 
planets in the mass range 1.65 Mj < M p sin(i) < 10 Mj is plotted 
for two intervals in semi-major axis. The solid histogram shows the 
distribution for planets with 0.1 AU < a < 1.185 AU, while the 
dashed histogram represents the distribution for 1.185 AU < a < 
2.5 AU. These distributions are not statistically distinguishable. 



mediately outside the gap for low mass planets that for 
high mass planets. If less mass piles up outside the tidal 
barrier, the torque will be reduced, and the migration 
velocity will go down. An example model of this kind, 
in which we assume that migration is fully suppressed 
for 0.5Mj planets but only partially suppressed for 2Mj 
planets, is also plotted in Figured] As expected (given 
the substantial difference between the relevant curves in 
Figure [5]) , this model displays order unity variations in 
the predicted slope of the mass function with radius. 

Observationally, an apparent paucity of high mass 
planets amongst the hot Jupiters has been noted by 
several authors (jZucker fc MazelJ [2002) . These plan- 
ets have presumably undergone interactions with their 
host stars (possibly including tidal circularization, mass 
loss, or stopping of migration due to entry into the stel- 
lar magnetosphere), any one of which could result in 
either mass loss or a mass-dependence to the planet's 
survival probability. At larger radii, however, there is 
no st rong evidence for any variation in the mass func- 
tion (|Marcv et al.l 12005). To illustrate this with current 
data, we c onsider the entire sample of known extraso- 
lar planets (|Butler et al.ll2006l ) with masses in the range 
1.65A/j < MpSin(i) < lOAfj, and orbital radii between 
0.1 AU < a < 2.5 AU (similar conclusions follow from 
analyses of smaller, more strictly selected samples). We 
divide this sample into 'inner' and 'outer' subsets, with 
equal numbers of planets in each. The dividing line be- 
tween the subsets falls at a semi- major axis of 1.185 AU. 
Figure [9] shows the cumulative distributions of M p sm(i) 
for these subsamples. By eye, and statistically, no sig- 
nificant differences are seen. This is consistent with the 
predictions of the simple migration model we have de- 
veloped here, and probably already rules out the large 
changes in the mass function with radius predicted in 
the case where high and low mass planets migrate in a 
qualitatively distinct manner. 
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Fig. 10. — Illustration of the effect of discontinuities in the 
efficiency of disk angular momentum transport on the predicted 
planet distribution. In this, not very realistic example, the viscosity 
is assumed to drop by a factor of 10 inside 1 AU (upper curve), 
leading to a large accumulation of planets at small orbital radii. As 
in Figure 6, the histogram shows the observed planet distribution 
and the lower curve the prediction of the fiducial migration model. 

5. SENSITIVITY TO STRUCTURE IN THE 
PROTOPLANETARY DISK 

The structure of the protoplanetary disk at radii 
of the order of 1 AU is uncertain theoretically, pri- 
marily as a consequence of the low ionization frac- 
tion which sows doubt as to the efficiency of angu- 
lar momentum t ransport driven by the magnetorota- 
tional instability dBalbus fe Hawlevlfl998U Gam:mIll996l: 
iTurner. Sano fe Dziourkevitchl 120*07 ). Moreover, most 
existing observatio nal constraint s are either explicitly 
(IWilner et al.ll2000ft or implicitly (jHartmann et al.lll998l : 
lArmitage. Clarke fe Pallal l2003i ) sensitive only to large 
radii in the disk. Not only is the slope of the steady-state 
surface density profile at small radii (E oc r~ 7 ) poorly 
known, but there could in principle be discontinuities in 
£ due to abrupt changes in the angular momentum trans- 
port efficiency (|Gammielll996l) . We have looked at how 
both of these possibilities affect the resulting distribution 
of planets. 

Holding other parameters of the model (primarily the 
planet formation radius) fixed, the resulting differential 
distributions of massive planets for different 7 are plot- 
ted in Figure O Here we have assumed a smooth surface 
density profile, which is fully characterized by the value 
of 7 which we vary between 7 = 0.5 and 7 = 1.5. The dif- 
ferences between these models at sub-AU radii are quite 
significant. If we consider the relative number of planets 
at 1 AU and 0.1 AU, the ratio varies from 2.4 (7 = 1.5) to 
3.8 (7 = 1) to 6.6 (7 = 0.5). The unbiased subsample de- 
fined earlier includes only a handful of planets at sub-AU 
radii, so comparison with existing data is not possible, 
but these differences are large enough that only a small 
unbiased sample of planets at these small radii would be 
needed to see variations of this magnitude. Changing 
the assumed planet formation radius has only a small 
influence on the distribution at sub-AU scales, so the in- 
fluence of different 7 is in principle separable from other 
unknown parameters of the model. 



Abrupt changes in the surface density profile — for ex- 
ample due to opacity transitions within the disk or rapid 
changes in the efficiency of angular momentum transport 
with radius, leave an even more distinctive fingerprint in 
the radial distribution of planets, provided only that mi- 
gration is not fully suppressed, so that r p is not linear in 
B. Such changes lead to a jump in the differential distri- 
bution of planets at the radius where the discontinuity 
occurs, and correspond to a change of slope in the cu- 
mulative distribution. As an example, Figure [TU] shows 
the predicted cumulative distribution if the viscosity is 
assumed to drop by a factor of 10 within 1 AU. Only 
rather unrealistic toy models of this kind (in which the 
jump is both large and situated squarely in the middle of 
the accessible radial range) can be ruled out using exist- 
ing data, but larger data sets could potentially constrain 
discontinuities in the disk physics quite well. 

6. CONCLUSIONS 

In this paper we have investigated the predictions of 
a rather simple, almost entirely analytic, model of giant 
planet formation and migration for the radial distribu- 
tion of massive extrasolar planets. Our main results are: 

1. The distribution of massive planets with radius, in- 
side the snow line beyond which planets are as- 
sumed to form, preserves information about the 
typical location of planet formation and the struc- 
ture of the protoplanetary disk at small radii. Sta- 
tistical analyses of large, uniformly selected sam- 
ples of massive planets could therefore in principle 
be used to recover information about planet forma- 
tion and disk physics. 

2. Abrupt changes in the disk properties with ra- 
dius, which may occur at small orbital radii, leave 
the most distinctive signature in the resulting 
planet distribution (a corresponding discontinuity 
in dN p /d\oga). Constraints on smooth surface 
density profiles at AU scales are also possible, but 
are likely to harder to obtain and more ambiguous. 

3. The mass function of extrasolar planets is not pre- 
dicted to vary detectably with radius as a conse- 
quence of migration, at least in the simplest mod- 
els. 

4. Existing data remains compatible with a sur- 
prisingly simple mo del, proposed previously 
( Armitag e et al.l [2002) , in which massive planets 
form beyond the snow line at a constant rate, and 
migrate inward through a smooth protoplanetary 
disk before becoming stranded when the gas disk 
is dispersed due to photoevaporation. However, the 
sample of planets suitable for a statistical compar- 
ison with the theory is limited, and as a result a 
wide class of alternative models remain viable or 
can only be ruled out with limited significance. 

Without a doubt, the simple theoretical model we have 
explored here is too simple. On a technical level, it would 
be desirable to extend the model to treat properly the 
case when planets form within a disk whose mass is only 
a small multiple of the planet mass, as this situation is a 
common outcome of the scenarios we favor. The model 
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also ignores some effects, such as mass flow across the 
gap and accretion onto the planet, that are almost cer- 
tainly important (ILubow. Seibert fc Artvmowicd 119991 : 
ILubow fc D'Aneelol l200(l ~ along with others, such 
as planet-planet scattering which has the side-effect 
of altering the orbital element s of surviving planets 
(|Ford. Havlickova fc Rasiol 12001). that could be signifi- 
cant. As such, the current agreement between our pre- 
dictions and the observed distribution of extrasolar plan- 
ets speaks more to the paucity of the data than to the 
validity of the theoretical model. Taking a broader view, 
however, the effects that we have ignored do not appear 
to be so intractable that they could not, in the future, be 
quantified and incorporated into a theoretical model of 
migration. Our results therefore imply that comparison 



of theoretical models to the larger samples of extrasolar 
planets that are forthcoming holds substantial promise 
for learning new information, about both planet forma- 
tion and protoplanetary disk physics, that is currently 
unavailable via other methods. 
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